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In  this  study,  we  are  simulating  the  turbulent  combustion  of  a  mixed  bluff-body  swirl  stabilized  flame  in 
a  gas  turbine  combustion  chamber  and  investigating  the  effects  of  different  parameters,  including  the 
swirl  number,  distance  between  the  air  and  fuel  nozzle  which  is  called  bluff  size,  equivalence  ratio,  inlet 
fuel  flow  rate,  and  the  inlet  air  velocity,  on  the  entropy  generation.  We  perform  the  process  of  the  design 
of  the  combustion  chamber  by  proposing  the  optimal  value  of  each  parameter  based  on  the  EGM  (en¬ 
tropy  generation  minimization)  method  under  the  two  maximum  allowable  temperature  and  size 
constraints.  Two  common  methods  of  entropy  generation  calculation,  one  based  on  the  overall  entropy 
balance  on  a  system  and  the  other  based  on  the  local  entropy  generation  rate  calculation,  are  used  and 
compared  in  this  study.  Our  results  show  that  the  deviation  between  the  total  entropy  generations 
calculated  by  the  two  methods  is  6.4%  in  average  which  is  an  acceptable  error  in  turbulent  combustion 
simulations.  Also,  the  two  opposing  factors,  namely  chemical  reaction  and  heat  transfer,  have  the  main 
contribution  to  the  total  entropy  generation. 

©  2015  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

The  efficient  utilization  of  energy  is  a  critical  concern  in  today's 
design  of  energy  systems.  Scientists  are  trying  to  approach 
maximum  efficiency  in  various  industrial  applications  as  much  as 
possible.  In  other  words,  it  is  favorable  to  have  minimum  irre¬ 
versibility.  In  order  to  achieve  this  goal,  the  exergy  analysis  and 
EGM  (entropy  generation  minimization)  method  1  can  be  applied 
to  optimize  the  design  of  energy  conversion  systems.  In  this 
method,  the  design  parameters  of  the  system  are  changed  and 
entropy  generation  is  computed  in  each  case.  The  optimum  design 
corresponds  to  the  case  with  the  minimum  entropy  generation 
while  some  additional  geometry  and  economic  constraints  should 
be  satisfied.  For  the  details  of  the  method,  one  can  consult  the  book 
by  Bejan  [1  .  Bejan  2]  also  reported  an  extensive  review  on  the 
applications  of  the  EGM  method  in  several  thermal  problems. 

The  EGM  method  is  vastly  incorporated  by  researchers  in 
different  applications.  Li  and  Faghri  [3]  performed  both  the  exergy 
balance  analysis  and  the  local  entropy  generation  analysis  to 
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investigate  the  irreversibility  within  fuel  cells.  Mahian  et  al.  [4] 
applied  EGM  to  study  the  effect  of  using  nanofluids  on  the  en¬ 
tropy  generation  between  two  concentric  rotating  cylinders  at 
constant  heat  flux  boundary  conditions.  They  found  that  the  en¬ 
tropy  generation  decreases  with  the  increase  of  volume  fraction  of 
nanoparticles.  Makhanlall  et  al.  [5  found  an  optimum  optical 
thickness  and  an  optimum  tilt  angle  where  the  heat  losses  are 
minimized  for  gas  filled  solar  collectors.  Escandon  et  al.  [6]  deter¬ 
mined  the  local  and  average  entropy  generation  rates  for  the  electro 
osmotic  flow  of  a  non-Newtonian  fluid  in  a  parallel  flat  plate  micro 
channel. 

EGM  method  has  also  been  implemented  in  several  heat 
transfer  studies  in  pipes.  For  example,  Ko  and  Ting  [7]  and  Ko  [8] 
reported  different  optimal  parameters  for  fully  developed  flow  in 
helical  pipes  at  the  constant  heat  flux.  Jarungthammachote  [9] 
investigated  entropy  generation  for  the  pipes  with  different 
shapes  of  cross  section,  and  Amani  and  Nobari  [10,11]  studied  en¬ 
tropy  generation  in  the  entrance  region  of  curved  pipes  both 
numerically  and  analytically. 

EGM  is  extensively  used  for  IC  (internal  combustion)  engine 
problems.  Caton  [12]  observed  that,  in  general,  availability 
destruction  decreases  when  combustion  temperature  is  increased. 
He  [13]  also  examined  the  availability  destruction  in  an  adiabatic, 
constant  volume  combustion  in  a  vessel.  He  discussed  the  effect  of 
combustion  temperature  and  equivalence  ratio  on  the  entropy 
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Nomenclature 

Z 

mixture  fraction 

2/ 

mass  fraction  of  ith  atomic  species 

X 

axial  coordinate 

Latin  symbols 

Xi 

mole  fraction  of  ith  species 

A 

area 

Cp 

specific  heat 

Greek  symbols 

D 

mass  diffusion  coefficient 

a 

inverse  effective  Prandtl  number 

E 

energy 

? 

bluff-size  ratio 

Gk 

production  of  turbulent  kinetic  energy 

e 

turbulent  energy  dissipation  rate 

h 

enthalpy 

4) 

equivalence  ratio 

h 

diffusion  flux  of  ith  species 

A 

thermal  conductivity 

k 

turbulent  kinetic  energy 

P 

dynamic  viscosity 

Mi 

molecular  weight  of  ith  species 

Pi 

chemical  potential  of  ith  species 

P 

pressure 

P 

density 

Q 

volume  flow  rate 

T 

stress  tensor 

R 

universal  gas  constant 

Ri 

mass  rate  of  creation  of  ith  species  by  chemical 

Subscripts 

reaction 

a 

air 

r 

radial  coordinate 

b 

bluff  body 

S 

swirl  number,  modulus  of  mean  rate  of  strain,  global 

ch 

due  to  chemical  reaction 

entropy  generation  rate 

eff 

effective 

Sij 

mean  strain  rate 

/ 

fuel 

Sc 

Schmidt  number 

gen 

entropy  generation  due  to  all  phenomena 

S'" 

local  entropy  generation  rate 

h 

due  to  heat  transfer 

s 

specific  entropy 

ij 

species  index,  coordinate  index 

T 

temperature 

m 

due  to  mass  transfer 

u 

velocity 

n 

normal  component 

Va 

mean  air-inlet  velocity 

t 

turbulent 

Yi 

mass  fraction  of  ith  species 

V 

due  to  viscous  dissipation 

generation.  Rakopoulos  and  Kyritsis  [14]  Compared  the  entropy 
generation  of  the  internal  combustion  engine  operation  for 
methane,  methanol,  and  dodecane  fuels.  Nakonieczny  [15  inves¬ 
tigated  entropy  generation  in  the  charge  exchange  system  of  a 
turbocharged,  intercooled  diesel  engine  and  discussed  the  effect  of 
various  variables,  such  as  the  intercooler  air  temperature,  air  pipe 
length,  and  timings  of  inlet  valve  closure  and  exhaust  valve  open¬ 
ing,  on  the  entropy  production.  Rakopoulos  and  Kyritsis  [16] 
studied  entropy  generation  during  the  combustion  of  the 
hydrogen-enriched  natural  and  landfill  gas.  Finally,  Rakopoulos  and 
Giakoumis  17]  presented  comprehensive  reviews  of  the  studies  on 
the  availability  analyses  of  IC  engines  published  up  to  2006.  They 
reported  fruitful  conclusions  about  the  availability  destruction, 
heat  transfer  rates,  and  availability  expelled  with  the  exhaust  gas. 
Also,  they  concluded  that  less  availability  is  destroyed  with  lighter 
fuels  such  as  methane  and  methanol  compared  to  the  longer-chain 
hydrocarbons  which  is  attributed  to  the  lower  entropy  generation 
of  the  decomposition  of  lighter  molecules.  Sieniutycz  [18]  investi¬ 
gated  the  relationship  between  entropy  generation  and  some  pa¬ 
rameters,  like  fuel  consumption  and  power  generation,  in  an 
engine.  The  availability  analysis  of  a  spark  ignition  bio¬ 
gas-hydrogen  engine  during  the  closed  part  of  the  engine  cycle 
was  performed  by  Rakopoulos  and  Michos  19].  Knizley  et  al.  [20] 
investigated  the  effects  of  the  fuel  type,  reactant  temperature, 
reactant  pressure,  fuel-air  equivalence  ratio  and  diluents  on  the 
entropy  generation  and  availability  destruction  in  a  constant  en¬ 
ergy-volume  (UV)  combustion  process. 

In  addition  to  IC  engines,  exergy  analysis  and  EGM  has  been 
applied  to  other  single-  and  multi-phase  combustion  processes. 
Dash  and  Som  [21  investigated  the  transport  processes  and  asso¬ 
ciated  irreversibility  of  droplet  combustion  in  a  convective  me¬ 
dium.  In  a  later  study,  Dash  and  Som  [22]  analyzed  thermodynamic 


irreversibility  of  the  spray  combustion  process  for  different  inlet 
pressures,  temperatures,  and  swirls.  Datta  [23]  studied  the  entropy 
generation  of  a  confined  laminar  diffusion  flame  and  the  influences 
of  the  inlet  air  temperature,  heat  loss  through  the  wall,  and 
equivalence  ratio  on  the  total  rate  of  entropy  generation.  Also, 
Nishida  et  al.  [24]  investigated  the  effects  of  these  parameters  on 
the  entropy  generation  of  both  premixed  and  diffusion  flames. 
Datta  [25]  extended  his  previous  study  [23]  to  account  for  the  role 
of  gravity  on  the  entropy  generation.  He  concluded  that  the  entropy 
generation  increases  by  decreasing  gravity.  Stanciu  et  al.  [26] 
investigated  the  irreversibility  in  both  laminar  and  turbulent 
diffusion  flames.  Two  combustion  models  are  used  by  them;  the 
multi-species  approach  based  on  the  eddy-break-up  model  for  the 
mean  reaction  rate,  and  the  assumed  probability  density  function 
for  a  conserved  scalar  that  relies  on  the  flame  sheet  model.  Yapici 
et  al.  [27-29]  investigated  the  effects  of  oxygen  fraction  in  the  air, 
equivalence  ratio,  and  swirl  on  the  entropy  generation  in  methane- 
and  hydrogen-fueled  combustion  chambers.  They  also  compared 
the  local  entropy  generation  for  different  fuels  [30].  Entropy  gen¬ 
eration  during  burning  of  a  spherical  fuel  droplet  was  calculated  by 
Raghavan  et  al.  31  ].  Investigation  of  the  influence  of  swirl  angle  on 
the  irreversibility  in  a  turbulent  diffusion  flame  was  performed  by 
Stanciu  et  al.  [32  .  Som  and  Datta  [33]  performed  a  fundamental 
study  on  the  thermodynamic  irreversibility  in  the  processes  of 
combustion  of  gaseous,  liquid,  and  solid  fuels.  They  discussed  the 
effects  of  the  free  stream  velocity,  particle  diameter,  ambient 
temperature,  and  gravity  on  the  entropy  generation  rate  in  both 
upward  and  downward  flows.  They  reported  the  decrease  of  the 
entropy  generation  by  reducing  the  Froude  number.  Analysis  of  the 
entropy  generation  of  laminar  propagating  hydrogen-enriched  co¬ 
flow  methane-air  triple  flames  was  conducted  by  Briones  et  al.  [34]. 
Chen  35  investigated  the  entropy  generation  of  laminar  premixed 
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hydrogen-air  combustion  of  planar  opposing  jets.  Chen  et  al.  [36] 
investigated  the  effects  of  hydrogen  addition  to  methane  on  the 
entropy  generation  in  a  similar  flame.  The  effect  of  equivalence 
ratio  and  hydrogen  addition  to  methane-air  laminar  counter-flow 
diffusion  flames  were  also  investigated  by  Chen  et  al.  [37,38  .  Bidi 
et  al.  [39]  studied  the  entropy  generation  of  combustion  in  porous 
media.  Sheikhi  et  al.  [40,41  ]  analyzed  local  entropy  generation  of 
turbulent  flows  using  LES  (large  eddy  simulation).  Emadi  and 
Emami  [42]  evaluated  entropy  generation  rate  and  exergy  loss  in  a 
turbulent  non-premixed  hydrogen-enriched  bluff-body  flame. 
Makhanlall  et  al.  [43]  performed  the  second  law  analysis  of  tur¬ 
bulent  diffusion  natural  gas  flames  and  investigated  the  importance 
of  radiative  energy  devaluation.  Entropy  generation  in  turbulent 
premixed  flames  was  analyzed  by  Farran  and  Chakraborty  [44] 
based  on  a  simple  chemistry  Direct  Numerical  Simulation  data¬ 
base  with  a  range  of  different  values  of  heat  release  parameter  and 
global  Lewis  number,  spanning  both  the  corrugated  flamelet  and 
thin  reaction  zone  regimes  of  combustion.  In  addition,  numerical 
simulation  of  biogas  flameless  combustion  was  performed  by 
Hosseini  et  al.  45]  to  analyze  the  entropy  generation  and  emission 
in  this  regime.  A  numerical  model  was  developed  for  the  entropy 
generation  and  exergy  destruction  analysis  of  hydrogen/air  pre¬ 
mixed  flames  in  micro-combustors  with  baffles  for  different  baffles 
height  and  mass  flow  rates  by  Jiang  et  al.  46].  Rana  et  al.  [47] 
evaluated  the  exergy  transfer  and  its  destruction  due  to  thermo¬ 
dynamic  irreversibility  in  the  process  of  flame  stabilization  in  a  heat 
recirculating  micro  combustor. 

In  the  present  study,  we  are  investigating  the  effects  of  different 
design  parameters  on  the  entropy  generation  of  methane-air 
combustion  in  a  gas  turbine  combustion  chamber.  This  involves 
the  study  of  two  important  geometrical  parameters,  i.e.  the  area  of 
the  air-inlet  nozzle  and  the  distance  between  the  air  and  fuel  nozzle 
called  bluff  size,  which  affect  the  fuel-air  mixing  process.  In  addi¬ 
tion,  the  impact  of  the  inlet  air  swirl  is  investigated.  Experiments 
show  that  these  parameters  have  considerable  influence  on  the 
mixing  and  pollutant  formation.  We  also  investigate  the  effect  of 
equivalence  ratio  and  its  interaction  with  the  other  parameters  on 
the  entropy  generation  rate.  The  optimal  design  is  proposed  based 
on  the  EGM  method  and  we  discuss  how  different  parameters 
affect  the  maximum  temperature  and  flame  height.  Two  common 
methods  of  entropy  generation  calculation,  one  based  on  applying 
the  entropy  balance  on  the  system  and  the  other  based  on  the  local 
entropy  generation  rate  calculation,  are  used  and  compared  in  this 
study.  Here,  we  call  the  first  method  the  direct  method  because  it  is 
based  on  direct  implementation  of  the  second  law  of  thermody¬ 
namics  on  a  system.  The  second  method  which  is  used  mostly  in 
CFD  (computational  fluid  dynamics)  studies  is  called  the  indirect 
method  here.  This  method  makes  it  possible  to  study  the  entropy 
generation  due  to  each  phenomenon  present  in  a  problem.  Yapici 
et  al.  [28]  performed  a  study  to  obtain  the  optimum  equivalence 
ratio  and  swirl  number,  but  he  only  considered  the  entropy  gen¬ 
eration  due  to  the  viscous  friction  and  heat  transfer,  while  we 
consider  the  entropy  generation  due  to  chemical  reaction,  mass 
transfer,  and  coupling  between  mass  and  heat  transfer  processes  in 
addition  to  those  factors  considered  by  Yapici  et  al.  28].  We  show 
that  not  only  are  these  entropy  generation  rates  important  in 
combustion  problems,  but  also  they  can  have  dominant  contribu¬ 
tion  to  the  total  entropy  generation  rate  and  cannot  be  neglected. 


combustion.  The  choice  of  these  models  is  based  on  the  trade-off 
between  accuracy  and  computational  cost  in  our  study  which  re¬ 
quires  the  simulation  of  a  large  number  of  combustion  chamber 
topologies  and  operating  conditions.  Although,  the  large  eddy 
simulation  offers  a  more  accurate  turbulence  closure,  its  usage  is 
limited  due  to  the  high  computational  costs.  Hence,  the  RANS 
(Reynolds-averaged  Navier-Stokes)  based  closures  have  still  the 
main  roll  in  practical  engineering  combustion  chamber  simulation 
and  design  48  .  For  the  problem  of  our  interest,  the  RNG  k-e  model 
with  the  modifications  accounting  for  the  swirl  is  a  wise  compro¬ 
mise  between  the  accuracy  and  computational  cost.  Here,  the  two- 
step  Westbrook-Drier  reaction  mechanism  [49]  which  is  appro¬ 
priate  for  fuel  lean  mixtures  governs  the  chemistry. 

In  the  following  paragraphs,  the  governing  equations  used  in 
our  study  are  introduced.  All  equations  are  given  using  Cartesien 
coordinate  notation  where  the  summation  convention  over 
repeated  indices  is  implied.  The  governing  equations  are  the 
transport  of  chemical  species  as 

s-yw  ■  -s* + *  (i» 

where  Uj  is  the  jth  component  of  the  (mean)  velocity  vector,  p  is  the 
density  of  the  mixture,  Y;  is  the  local  mass  fraction  of  species  i,  R\  is 
the  mass  rate  of  creation  or  depletion  of  species  i  by  the  chemical 
reaction,  and  Ji  is  the  total,  molecular  and  turbulent,  flux  of  species  i 
and  is  modeled  by 

av. 

Jij  =  - PDeff ^  (2) 

where  Deff  is  the  effective  mass  diffusivity  of  each  species  in  the 
mixture  and  defined  in  the  RNG  model  as 


pa 


and  a  is  computed  by  the  following  equation 


a  —  \  .3929 

0.6321 

a  +  2.3929 

0.3679 

 P 

ao  —  1 .3929 

a$  ~T  2.3929 

Mf 

(3) 


(4) 


where  ao  =  1 1  Sc  and  Sc  is  the  molecular  Schmidt  number  defined  by 

Sc^  (5) 

pD 


and  D,  p,  p,  and  pieff  are  the  molecular  mass  diffusivity,  density, 
molecular  viscosity,  and  effective  viscosity  of  the  mixture, 
respectively. 

The  momentum  equation  is  written  by 


dP  0Ty 

dxf  ^  dXj 


where  P  is  the  pressure  and  Ty  is  the  effective  stress  tensor. 


ry  —  ^ eff 


2  dU[ 

3 


2.  Mathematical  model 


The  energy  conservation  equation  is 


2.1.  Governing  equations 

In  this  study,  we  use  RNG  (renormalization  group)  k-e  model  for 
the  turbulence  closure  and  the  eddy-dissipation  model  for  the 


0 

dXj 


Uj(pE  +  P)]  -  ^  (^effQX_  ^  hJi.j  +  UjTij^J  + 


(8) 


where  Sh  =  -  includes  the  heat  of  chemical  reaction, 
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r  P  UJUj 

E  =  h - h^r1 

P  2 


(9) 


h  is  the  sensible  enthalpy,  and  hjh  the  formation  enthalpy  of  species 
i.  Xeff  is  the  effective  conductivity  and  is  given  by 


A eff  =  aCpPeff  (10) 

where  a  in  Eq.  (10)  is  computed  from  Eq.  (4)  substituting  ao  =  A//xCp. 
A  is  the  molecular  conductivity  and  Cp  is  the  mixture  specific  heat 
capacity. 

Two  additional  equations  in  RNG  k— e  turbulence  model  are 


0 

Wj 


d  (  dk\ 

=w][a^ffwj) 


(ii) 


4 {pUje) = i  {aefleffSj) +  i (CleCfc ' Clepe)  ~x  ( 1 2) 

where  k  is  the  turbulence  kinetic  energy  and  e  is  its  dissipation  rate. 
ak  and  a£  are  the  inverse  effective  Prandtl  numbers  for  k  and  e, 
respectively,  and  are  calculated  from  Eq.  (4)  substituting  ao  =  1.  The 
effective  viscosity  in  RNG  k-e  turbulence  model  is  computed  by 


Me//  =  M  +  Mt, 


Mt  pQi 


and 


(13) 


(14) 


X  =  Epp 


1  +  /??73  k 


(15) 


r)  =  S—  (16) 

£ 

where  Sy  is  the  mean  rate  of  strain  tensor  and  defined  as 


_  1  (dUi  duj\ 
ij  ~  2  (ox,  +  dXi) 

and  its  modulus,  S,  is  given  by 


(17) 


Table  1 

The  molecular  properties  of  species  used  in  our  study.  D,  =  2.88  x  10-5m2/s  for  all 
species. 


Species 

Pi  (kg/ms) 

A*  (w/m  K) 

ch4 

0.0332 

1.087  x  10~5 

o2 

0.0246 

1.919  x  10  5 

co2 

0.0145 

1.370  x  10^5 

CO 

0.0250 

1.750  x  10~5 

h2o 

0.0261 

1.340  x  10  5 

n2 

0.0242 

1.663  x  10  5 

smaller  than  their  turbulent  counterparts.  The  sensible  enthalpy  is 

‘T 


related  to  the  mixture  specific  heat  via  h 


CpdT. 


! ref 


2.2.  Entropy  generation  calculation 

Here,  we  calculate  the  entropy  generation  using  two  approaches. 
In  the  first  method,  which  is  called  “direct  method”  hereafter,  we 
use  the  second  law  directly  to  compute  the  total  entropy  generation 
rate  in  a  system,  i.e.  for  an  adiabatic  combustion  chamber  per¬ 
forming  no  work  on  its  surrounding;  the  total  rate  of  entropy 
generation  equals  the  net  outflow  of  entropy  from  the  system 

Y)  y,s,j  u„dA 
(19) 

Therefore,  Sgen  is  computed  by  integrating  mass  flux  of  entropy 
over  the  inlet  and  outlet  boundaries.  Note  that,  s/  is  the  entropy  of 
species  i  in  the  mixture,  i.e. 

sf  =  Sj(7\  Pi)  =  s?(T)  -  ^ In  00  (20) 

where  R  is  the  universal  gas  constant. 

In  the  second  method,  called  “indirect  method”,  we  calculate 
the  local  entropy  generation  rates  (per  unit  volume)  due  to 
different  phenomena  present  in  our  problem.  The  causes  of  entropy 
generation  in  a  single-phase  combustion  problem  are  the  viscous 
friction,  heat  transfer,  mass  diffusion,  coupling  term  between  the 
mass  and  heat  transfer,  and  chemical  reaction  [52  .  These  rates  can 
be  written  as 


■>gen 


Sout  —  Sin  = 


undA-  A 


in 


rfff  _  T_g_  ( dtlj\ 

v  ~  T  VaxJ 


(21a) 


S=V/2S^  (18) 

The  constants  of  the  model  are  CM  =  0.0845,  C\e  =  1.42, 
C2e  =  1.68, 7]o  =  4.38  and  (3  =  0.012. 

The  mixture  density  is  calculated  using  the  ideal  gas  law. 
Combining  the  properties  of  different  species,  the  properties  of  the 
mixture  are  calculated;  for  the  mixture  specific  heat,  the  linear  rule, 

i.e.  Cp  =  YiCp,u  is  used  where  Cp/s,  as  well  as  sf's  in  Eq.  (20),  of 

i 

different  species  are  determined  as  a  function  of  temperature  from 
the  7-coefficent  NASA  polynomials  [50];  and  for  the  transport  co¬ 
efficients  of  the  mixture,  including  molecular  viscosity,  diffusivity, 
and  conductivity,  the  “Wilke  rule”  [51  is  used  where  the  molecular 
transport  properties  of  each  species  are  assumed  constant  and 
given  in  Table  1.  Note  that  the  constant  molecular  transport  prop¬ 
erties  assumption  has  a  negligible  effect  on  the  results  since  these 
properties  are  about  two,  or  even  more,  orders  of  magnitude 


<jw  =  hn( EL  EL) 

h  T2  {dXjdXjJ 


(21b) 


<-///  _  l,RDeff  ( 9 Yi  (Tf.) 

m  “  4^  YiMi  \dxj  dx)J 


(21C) 


(21  d) 


s'"h  =  -Y^r  (21e) 

i 

where  m  is  the  chemical  potential  of  species  i  in  the  mixture  which 
is  given  by 
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Hi  =  h°f  i  +  h,  -  Tst  (22) 

The  rates  in  Eq.  (21a— e)  are  also  given  in  axisymmetric  co¬ 
ordinates  in  Appendix  A.  The  total  entropy  generation  rate  per  unit 
volume,  Sgen,  is  the  sum  of  all  local  rates  in  Eq.  (21).  Eventually,  the 
total  entropy  generation  rate  in  the  system  is  obtained  by  inte¬ 
grating  S'"  over  the  entire  system. 

Sgen  =  j  S’"endV  (23) 

V 


3.  Numerical  method 

The  geometry  of  the  combustion  chamber  and  a  typical  grid 
used  in  this  study  are  shown  in  Fig.  1.  r0  is  the  radius  of  the  com¬ 
bustion  chamber,  r\  the  outer  radius  of  the  air-inlet  nozzle,  X2  the 
inner  radius  of  this  nozzle,  and  r 3  the  radius  of  the  fuel  inlet  nozzle. 

For  the  numerical  discretization  of  the  governing  equations, 
finite  volume  method  is  applied.  The  diffusion  terms  are  discretized 
by  the  second-order  central  scheme  and  convection  terms  by  the 
standard  second-order  upwind  scheme  [53].  The  SIMPLE  algorithm 
[54]  is  used  for  the  pressure-velocity  decoupling  of  the  Navier — 
Stokes  equations.  The  standard  non-linear  limiter,  which  is  derived 
from  the  work  of  Barth  and  Jespersen  [55],  is  incorporated  to  pre¬ 
serve  the  monotonicity  of  the  convection  discretization  scheme. 

3.1.  Verification  and  validation 

To  report  the  grid  sensitivity  test,  the  profile  of  different  pa¬ 
rameters  in  various  cross  sections  of  the  chamber  are  depicted  in 
Figs.  2-4.  In  these  figures,  one  case  corresponds  to  the  typical  grid 
shown  in  Fig.  1  and  the  other  case  corresponds  to  a  refined  grid 
obtained  by  dividing  each  computational  cell  of  the  typical  grid  to 
four  cells.  It  is  observed  that  the  results  of  different  grids  are  in  fine 
agreement  with  each  other  and,  therefore,  the  grid  independency 
of  the  solution  is  confirmed. 

In  order  to  check  our  entropy  generation  calculations,  we  also 
verify  our  numerical  computation  for  a  special  case,  i.e.  the  laminar 
flow  in  a  straight  duct  with  constant  heat  flux  at  the  wall  which  has 
an  analytical  solution  for  the  profiles  of  the  entropy  generation  rate 
due  to  the  viscous  friction  and  heat  transfer  at  the  fully  developed 
cross  section.  In  this  case,  no  estimation  is  required  to  compute  the 
entropy  generation  in  contrast  to  the  turbulent  flame  cases 
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Fig.  1.  a)  The  geometry  of  combustion  chamber  and  b)  a  part  of  combustion  chamber 
near  the  inlet  and  the  typical  (146  x  400)  grid  used  for  the  simulations. 


reported  in  our  study.  Therefore,  the  only  source  of  error  is  nu¬ 
merical  errors  which  can  be  evaluated  through  this  verification.  The 
comparison  is  illustrated  in  Fig.  5  for  water  flow  in  a  duct  with 
radius  ro  =  6  mm,  ito  =  0.01  m/s,  To  =  273  K,  and  q  =  1000W/m2. 
The  analytical  solutions  for  the  entropy  generation  rate  profiles  are 


rfff  _  H 

h  ~  AT2 


2a 

fimX 0. 


1  2 


+ 


(24) 


(25) 

where  a  =  A/pcp  and  Tm  and  um  are  the  bulk  temperature  and  ve¬ 
locity  of  the  flow  in  each  cross  section.  As  it  can  be  observed  in 
these  figures,  our  calculations  are  in  close  agreement  with  the 
analytical  solutions. 

Subsequently,  in  order  to  validate  the  present  model,  we 
simulate  the  well-known  methane-air  turbulent  non-premixed  jet 
flame  “Sandia  flame  D”  and  compare  the  results  of  our  model  with 
the  experimental  measurements  56].  This  flame  has  also  been 
used  for  entropy  generation  validation  by  Sheikhi  et  al.  [40  .  For  the 
details  of  the  burner  geometry  and  inlet  boundary  profiles,  readers 
can  refer  to  the  corresponding  web  page  [56  .  We  use  the  same 
boundary  conditions  and  computational  grid  as  our  previous  study 
[57]  while  the  present  model  differs.  The  mean  entropy  is  extracted 
from  the  measured  instantaneous  temperature  and  species  mass 
fractions  [40,56  .  As  it  is  observed  in  Fig.  6,  the  computed  results 
have  in  average  9.5%  overprediction  with  respect  to  the  experiment 
which  is  within  the  acceptable  range  of  the  errors  of  estimations 
described  in  the  following  section.  This  deviation  in  entropy  pre¬ 
diction  is  mainly  due  to  temperature  overprediction,  typical  of 
RANS/Eddy-dissipation  models  [57  .  Note  also  that  neglecting  the 
radiation  effect  generally  results  in  an  overprediction  of  the  flame 
temperature,  but  for  such  a  non-sooty  methane-air  flame,  i.e. 
Sandia  D,  this  simplification  does  not  lead  to  significant  errors 
(about  8%  overprediction  as  investigated  by  Xu  et  al.  [58  ). 


4.  Results  and  discussion 


4.2.  Case  studies 


The  bluff  size,  or  the  width  of  the  wall  between  the  air  and  fuel 
nozzle,  has  an  important  effect  on  the  mixing  and  efficiency  of  the 
combustion  in  a  gas  turbine  combustion  chamber.  To  study  the 
effect  of  the  bluff  size,  we  define  a  new  dimensionless  parameter  as 

Z  =  T-  (26) 


where  Aa  is  the  area  of  the  air-inlet  nozzle  and  At  is  the  area  of  the 
annular  wall  (bluff)  between  the  air-inlet  and  fuel-inlet  nozzles. 
Hereafter,  f  is  called  bluff-size  ratio.  Simulations  are  carried  out  for 
the  bluff-size  ratio  in  the  range  0.5-10.  The  other  important  factor 
is  the  air-inlet  swirl.  The  well-known  swirl  number  parameter,  S, 
which  is  the  ratio  of  tangential  to  axial  momenta 


■n 


uxuer2dr 


r2 


X\  1  u2xdx 


(27) 


is  used  to  quantify  the  swirl  intensity.  Here,  we  study  the  effect  of 
swirl  number  in  the  range  0-0.4.  The  3rd  design  parameter  is  the 
air-inlet  velocity,  Va,  which  can  be  changed  by  altering  the  air-inlet 
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(a)  (b)  (c) 

Fig.  2.  The  profiles  of  the  local  entropy  generation  rate  due  to  chemical  reaction  at  different  cross  sections  using  the  typical  grid  (146  x  700)  and  the  refined  grid  (292  x  800)  a) 
x  =  0.4  cm,  b)  x  =  1  cm,  and  c)  x  =  10  cm. 
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Fig.  3.  The  profiles  of  the  local  entropy  generation  rate  due  to  heat  transfer  at  different  cross  sections  using  the  typical  grid  (146  x  700)  and  the  refined  grid  (292  x  800)  a) 
x  =  0.4  cm,  b)  x  =  1  cm,  and  c)  x  =  10  cm. 
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Fig.  4.  The  profiles  of  the  temperature  at  different  cross  sections  using  the  typical  grid  (146  x  700)  and  the  refined  grid  (292  x  800)  a)  x  =  0.4  cm,  b)  x  =  1  cm,  and  c)  x  =  10  cm. 
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r/D 

Fig.  5.  The  local  entropy  generation  rate  at  section  x/D  =  65;  up:  due  to  heat  transfer, 
down:  due  to  viscous  friction. 


nozzle  area  assuming  constant  air  flow  rate.  We  consider  different 
values  of  the  air-inlet  velocity  in  the  range  3.969-31.752  m/s  in  our 
simulations.  The  other  important  design  parameter  is  the  equiva¬ 
lence  ratio,  </>,  which  indicates  the  ratio  of  air  to  fuel  mass  flow  rates. 
This  parameter  takes  different  values  in  the  range  0.5-1.4  in  our 
simulations.  We  also  consider  different  working  condition  of  the 
combustion  chamber  by  taking  into  account  the  different  values  of 
fuel-inlet  flow  rate  in  the  range  5-10  1pm  (liter  per  minute).  Other 
geometrical  parameters  of  the  combustion  chamber  (see  Fig.  1 )  are 
r0  =  6  cm,  r3  =  0.4  cm,  L  =  50  cm.  Note  that,  \\  and  r2  take  different 
values  in  the  range  0.714-1.721  cm  and  0.438-1.197  cm,  respec¬ 
tively,  for  introducing  different  ?  or  Va. 

4.2.  Entropy  generation  minimization 

As  it  was  mentioned  earlier,  using  indirect  method  enables  us  to 
study  the  effect  of  different  phenomena  on  the  entropy  generation 
separately.  Beside  this  fact,  using  these  two  methods  for  the  total 
entropy  generation  calculation  and  comparing  them  has  another 
advantage.  To  explain  this  advantage,  note  that  the  formulae  of 
both  direct  method,  Eq.  (19),  and  indirect  method,  Eq.  (21),  include 
the  multiplication  of  different  variables,  i.e.  non-linear  terms.  Using 
these  formulae  for  turbulent  flows  replacing  the  instantaneous 
variables  with  their  means  and  molecular  transport  coefficients,  fi, 
A,  and  D,  with  the  effective  transport  coefficients,  geff,  Aejy,  and  Dejj, ,  is 
an  approximation.  Therefore,  computing  and  comparing  the  total 


Fig.  6.  Comparison  of  the  present  model  prediction  and  experimental  data  (Ref.  [56  ). 
The  radial  profile  of  the  mixture  specific  entropy,  s  =  two  cross  sections;  x/ 

D  =  7.5  (up)  and  x/D  =  15  (down).  D  =  7.2  mm  is  the  fuel  nozzle  diameter.  Note  that  to 
determine  the  total  entropy  generation  rate  between  the  inlet  and  a  specific  cross 
section,  this  profile  is  integrated  using  Eq.  (19)  over  that  cross  section,  e.g.  Aout:  x/ 
D  =  7.5  or  A0ut:  x/D  =  15. 


entropy  generation  using  these  two  independent  approximations 
would  be  extremely  useful  to  validate  conclusions  extracted  from 
our  results.  Note  also  that  using  such  approximations  for  calcu¬ 
lating  the  entropy  generation  in  turbulent  flows  is  unavoidable. 
This  is  because  the  instantaneous  formulation,  see  Eq.  (21 ),  involves 
the  gradient  of  different  variables.  Therefore,  the  exact  calculation 
of  the  entropy  generation  requires  the  usage  of  two-point  PDF 
(probability  density  function)  methods  which  include  the  gradients 
of  all  variables,  in  addition  to  the  variables  themselves,  as  the  state 
variables  and  is  beyond  the  current  scope  of  turbulent  combustion 
modeling  and  simulation. 

Fig.  7  shows  the  total  entropy  generation  calculated  using  both 
direct  and  indirect  methods  for  quite  a  lot  of  simulation  cases 
involving  different  equivalence  ratios  and  inlet-fuel  flow  rates.  The 
results  of  both  methods  are  in  satisfactory  agreement  especially 
where  </>  <  0.7.  The  maximum  deviation  between  the  two  methods 
reaches  about  16%  at  the  fuel-lean  limit.  The  deviation  between  the 
two  methods  is  depicted  in  Fig.  8.  The  average  deviation  for  all 
simulations  is  6.4%  which  is  acceptable  in  turbulent  combustion 
modeling  and  supports  the  assumptions  underlying  the  entropy 
generation  calculations  in  our  study.  As  it  is  observed  in  Fig.  7,  the 
total  entropy  generation  reduces  by  increasing  the  equivalence 
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Fig.  7.  The  entropy  generation  rates  for  various  equivalence  ratios  and  inlet-fuel  flow 
rates.  Note  that  1pm  in  the  figures  refers  to  liter  per  minute.  The  other  parameters  are 
S  =  0,  Aa  =  2  cm2,  £  =  10,  and  Tiniet  =  298.15  K. 
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Fig.  9.  The  entropy  generation  rate  due  to  different  phenomena  versus  equivalence 
ratio.  The  other  parameters  are  Qf  =  10  1pm,  S  =  0,  Aa  =  2  cm2,  £  =  10,  and 
Tiniet  =  298.15  K. 


ratio  and  an  optimum  case  cannot  be  proposed  based  on  EGM.  This 
trend  was  also  reported  by  Yapici  et  al.  [28  in  spite  of  the  fact  that 
they  only  considered  viscous  and  heat  transfer  entropy  generation 
rates.  This  behavior  can  be  explained  via  Fig.  9.  This  figure  shows 
that  the  entropy  generations  due  to  chemical  reaction  and  heat 
transfer  have  the  dominant  contributions  to  the  total  entropy 
generation  in  the  combustion  chamber.  This  conclusion  is  also 
confirmed  by  the  next  results  and  diagrams  presented  in  this  paper. 
According  to  Fig.  9,  increasing  the  equivalence  ratio  towards  the 
fuel-rich  region  causes  the  heat  transfer  and  especially  chemical 
entropy  generations,  i.e.  the  dominant  parts  of  the  total  entropy 
generation,  to  decrease  monotonically.  This  is  due  to  the  drop  of 
combustion  intensity  and  the  resulted  heat  release  by  enriching  the 
fuel  in  the  combustion  chamber. 

Therefore,  the  optimal  value  of  the  equivalence  ratio  cannot  be 
chosen  only  based  on  minimizing  the  entropy  generation  and  is 
usually  selected  based  on  the  combustion  efficiency  and  the  MAT 
(maximum  allowable  temperature)  in  the  combustion  chamber. 
However,  the  MAT  criterion  and  other  criteria,  e.g.  geometrical, 
economical,  etc.,  can  be  regarded  as  the  constraints  in  the  EGM 
method  which  results  in  an  optimization  problem  under  several 
constraints.  Note  that,  the  equivalence  ratio  is  the  most  effective 
parameter  of  controlling  the  maximum  temperature  in  the  com¬ 
bustion  chamber.  This  can  also  be  seen  in  our  results  by  comparing 
Fig.  10  and  the  results  of  the  next  parts.  The  MAT  criterion  is 
important  from  two  different  aspects;  one  corresponds  to  the 
thermal  resistance  of  materials  used  to  fabricate  the  combustion 
chamber  and  other  downstream  parts,  e.g.  turbine  blades.  The 
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Fig.  8.  The  deviation  between  direct  and  indirect  methods  of  entropy  generation 
calculation. 


other  concern  is  related  to  the  pollutant  emission  control  because 
this  temperature  is  the  prominent  factor  affecting  the  NOX  (mono¬ 
nitrogen  oxides)  production.  To  account  for  the  MAT  constraint,  the 
maximum  temperature  within  the  combustion  chamber  is  also 
reported  in  Fig.  10.  For  example,  if  MAT  equals  1900  K,  the  value  of 
equivalence  ratio  is  limited  to  </>  =  0.5.  This  choice  also  corresponds 
to  a  fuel-lean  condition  typical  of  practical  combustion  chambers. 
Beside  MAT  constraint,  we  consider  another  constraint  regarding 
the  size  of  the  combustion  chamber  which  is  discussed  separately 
in  section  4.3.  We  also  checked  all  of  our  optimizations  presented  in 
this  part  to  satisfy  this  second  constraint.  Therefore,  our  optimi¬ 
zations  are  actually  performed  under  two  constraints. 

In  the  next  step,  we  investigate  the  effect  of  the  bluff-size  on  the 
entropy  generation  and  try  to  propose  a  superior  design  based  on 
EGM  method.  The  total  entropy  generation  within  the  combustion 
chamber  is  reported  for  different  bluff-size  ratios  in  Fig.  11  using 
both  direct  and  indirect  methods.  The  maximum  temperature 
within  the  combustion  chamber  is  also  shown  in  this  figure  to 
perform  an  optimization  under  the  MAT  constraint.  The  deviation 
between  the  direct  and  indirect  methods  is  about  13%  which  is 
acceptable.  Note  that  the  average  deviation  is  far  below  this  value 
because  for  this  case;  </>  =  0.5  which  is  in  fuel-lean  region  (refer  to 
Fig.  8).  Another  interesting  point  which  is  inferred  from  Fig.  11  is 
that  in  spite  of  some  deviation  between  the  total  entropy  genera¬ 
tions  computed  by  the  two  methods,  the  optimum  condition  pro¬ 
posed  by  both  of  them  is  approximately  identical.  In  other  words, 
both  methods  result  in  the  same  optimal  design.  This  further 
supports  the  optimal  conditions  proposed  in  our  study.  As  it  can  be 
observed  in  Fig.  11,  the  optimum  condition,  corresponding  to  the 


Fig.  10.  The  maximum  temperature  within  the  combustion  chamber  for  various 
equivalence  ratios  and  inlet-fuel  flow  rates.  The  other  parameters  are  S  =  0,  Aa  =  2  cm2, 
l  =  10,  and  Tiniet  =  298.15  K. 
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Fig.  11.  The  total  entropy  generation  versus  the  bluff-size  ratio  computed  using  direct 
and  indirect  methods  (up),  and  maximum  temperature  inside  the  combustion  cham¬ 
ber  versus  the  bluff-size  ratio  (down).  The  other  parameters  are  Qf  =  10  1pm,  S  =  0, 
Aa  =  2  cm2,  0  =  0.5,  Va  =  15.876  m/s,  and  Tiniet  =  298.15  K. 

minimum  total  entropy  generation  under  the  constraint  of 
MAT ~  1900/C,  is  at  ?  =  2.  The  next  figures  explain  this  result  further. 
Fig.  12  shows  the  entropy  generation  rate  due  to  different  phe¬ 
nomena  separately.  According  to  this  figure,  the  entropy  generation 
due  to  chemical  reaction  increases  by  the  increase  of  the  bluff-size 
ratio  while  the  heat  transfer  entropy  generation  decreases.  These 
two  phenomena  make  the  main  contribution  to  the  total  entropy 
generation.  This  behavior  is  explained  through  Figs.13-15.  For 
three  values  of  the  bluff-size  ratio,  the  flow  streamlines  are 
depicted  in  Fig.  13,  and  the  contours  of  the  fuel  mass  fraction  and 


0.0  2.0  4.0  6.0  8.0  10.0 

Fig.  12.  The  entropy  generation  rate  due  to  different  phenomena  versus  the  bluff-size 
ratio.  The  other  parameters  are  Qf  =  10  1pm,  S  =  0,  Aa  =  2  cm2, 0  =  0.5,  Va  =  15.876  m/s, 
and  Tiniet  =  298.15  K.  Line  legend:  the  same  as  Fig.  9. 
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temperature  are  depicted  in  Figs.  14  andl5,  respectively.  Note  that 
these  figures  show  a  part  of  the  combustion  chamber  near  the  in¬ 
lets.  For  example,  in  Fig.  13,  one  can  see  the  eddies  generated  in  the 
central  recirculating  region  of  the  chamber  near  the  inlet  streams. 
There  is  also  another  eddy  which  generated  near  the  corner  of  the 
combustion  chamber  due  to  the  flow  expansion  and  is  out  of  the 
depicted  frame  in  Figs.  13-15.  This  region  is  called  the  corner 
recirculating  region  in  a  combustion  chamber.  For  the  sake  of 
completeness,  this  corner  region  is  shown  in  Fig.  16  where  you  can 
see  the  corner  eddy.  The  eddies  in  central  recirculating  region  play 
an  important  role  in  fuel-air  mixing.  As  it  is  observed  in  Fig.  13,  by 
decreasing  the  bluff-size  ratio,  i.e.  increasing  the  bluff  area  with 
respect  to  the  air-inlet  area,  the  number  and  complexity  of  the 
generated  eddies  increases.  This  leads  to  the  increase  in  the  mixing 
rate  which  reduces  the  chemical  entropy  generation  rate  while  the 
entropy  generation  rate  due  to  heat  transfer  rises  due  to  the  in¬ 
crease  in  the  reaction  heat  release  (see  Fig.  12).  In  other  words,  the 
chemical  entropy  and  heat  transfer  entropy  are  opposing  phe¬ 
nomena  in  the  optimization  of  the  bluff-size  ratio.  The  trade-off 
between  these  opposing  rates  occurs  at  ?  =  2  corresponding  to 
the  minimum  total  entropy  generation.  Considering  the  maximum 
temperature  versus  the  bluff-size  ratio  in  Fig.  11,  this  optimum 
condition  also  satisfies  the  MAT  constraint. 


Fig.  13.  Streamlines  in  a  part  of  the  combustion  chamber,  near  the  inlets,  for  a)  £  =  0.5, 
b)  £  =  2  (optimized  case),  and  c)  £  =  10.  The  other  parameters  are  Qf=  10  1pm,  S  =  0, 
Aa  =  2  cm2,  0  =  0.5,  Va  =  15.876  m/s,  and  Txn/et  =  298.15  K. 
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Fig.  15.  The  temperature  distribution  in  the  combustion  chamber,  for  a)  £  =  0.5,  b) 
£  =  2  (optimized  case),  and  c)  £  =  10.  The  other  parameters  are  Qf  =  10  1pm,  S  =  0, 
Aa  =  2  cm2,  0  =  0.5,  Va  =  15.876  m/s,  and  Tiniet  =  298.15  K. 


(c) 

Fig.  14.  The  fuel  mass  fraction  in  a  part  of  the  combustion  chamber,  near  the  inlets,  for 
a)  £  =  0.5,  b)  £  =  2  (optimized  case),  and  c)  £  =  10.  The  other  parameters  are 
Qf  =  10  1pm,  S  =  0,  Aa  =  2  cm2,  0  =  0.5,  Va  =  15.876  m/s,  and  T/n;et  =  298.15  K. 


Next  stage  is  to  investigate  the  effect  of  swirling  flow  on  the 
performance  of  the  combustion  chamber  using  EGM  method.  Both 
total  entropy  generation  and  maximum  temperature  within  the 
combustion  chamber  for  various  swirl  numbers  are  depicted  in 
Fig.  17.  As  it  can  be  seen  in  this  figure,  the  total  entropy  generation 
increases  very  slightly  by  increasing  the  swirl  number.  In  fact,  by 
varying  the  swirl  number  in  the  range  of  0-0.4,  the  total  entropy 
generation  changes  by  only  4%  for  £  =  2  and  </>  =  0.5.  On  the  other 
hand,  the  maximum  temperature  increases  by  increasing  the  swirl 
number  violating  the  MAT  criterion.  Therefore,  the  proposed  value 
of  the  swirl  number  is  S  =  0  here.  Note  that,  due  to  the  limitation  of 
our  model,  only  the  low  to  moderate  swirl  numbers  are  considered 
here  which  combine  with  the  recirculating  flow  generated  by  the 
presence  of  the  bluff.  The  variation  of  the  entropy  generation  rate 
due  to  different  phenomena  by  the  change  of  the  swirl  number  is 
also  shown  in  Fig.  18.  Increasing  the  swirl  number  causes  the  in¬ 
crease  of  the  heat  transfer  entropy  and  decrease  of  chemical  en¬ 
tropy  generations.  This  variation  is  more  tangible  in  the  range 
S  <  02  and  after  that  it  is  not  much  recognizable. 

The  final  step  is  to  study  the  effect  of  air-inlet  velocity  on  the 
performance  of  the  combustion  chamber.  The  total  entropy  gen¬ 
eration  for  various  air-inlet  velocities  is  demonstrated  in  Fig.  19.  As 
it  is  observed  in  this  figure,  the  total  entropy  generation  reduces 


slightly  by  decreasing  the  air-inlet  velocity,  Va.  By  varying  this 
parameter  in  the  range  of  3.969-31.752  m/s,  the  total  entropy 
generation  changes  6%  for  f  =  2,  S  =  0  and  </>  =  0.5.  In  addition, 
Fig.  19  shows  that  the  maximum  temperature  in  the  chamber  is 
approximately  constant  when  Va  is  larger  than  15.876  m/s  and  in¬ 
creases  sharply  by  decreasing  its  value  beyond  15.876  m/s.  Thus,  in 
spite  of  the  (slightly)  smaller  value  of  entropy  generation  for 
smaller  air-inlet  velocities,  the  optimum  value  of  this  parameter  is 
Va  =  15.876  m/s  to  satisfy  the  MAT  constraint.  Fig.  20  shows  the 
entropy  generation  rate  due  to  each  factor  separately.  It  demon¬ 
strates  that  increasing  the  inlet  air  velocity  causes  the  increase  of 
the  heat  transfer  entropy  and  decrease  of  chemical  entropy 


Fig.  16.  Streamlines  in  the  corner  of  the  combustion  chamber  (corner  recirculating 
region)  corresponding  to  the  case  of  Fig.  13(b). 
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Fig.  17.  The  total  entropy  generation  versus  the  swirl  number  computed  using  direct 
and  indirect  methods  (up),  and  maximum  temperature  inside  the  combustion  cham¬ 
ber  versus  the  swirl  number  (down).  The  other  parameters  are  Qf  =  10  1pm,  £  =  2, 
Aa  =  2  cm2,  </>  =  0.5,  Va  =  15.876  m/s,  and  Tin;et  =  298.15  K. 

generations.  This  change  is  sharper  for  cases  with  lower  air-inlet 
velocities  than  15  m/s  and  after  that  it  is  not  much  tangible. 

4.3.  The  flame  height  as  the  second  constraint 


Fig.  19.  The  total  entropy  generation  versus  the  air-inlet  velocity  using  direct  and 
indirect  methods  (up),  and  maximum  temperature  inside  the  combustion  chamber 
versus  the  air-inlet  velocity  (down).  The  other  parameters  are  Qf  =  10  1pm,  £  =  2, 
Aa  =  2  cm2,  S  =  0,  </>  =  0.5,  and  T,-n/et  =  298.15  K. 

Here,  the  flame  height  is  calculated  based  on  the  mixture  frac¬ 
tion  definition  with  Bilger's  coefficients  as 


In  practical  design  of  energy  systems,  the  geometrical  con¬ 
straints  have  to  be  considered  almost  always.  One  of  the  important 
factors  which  determine  the  size  of  a  combustion  chamber  is  the 
flame  height  because  it  is  desirable  that  the  combustion  process 
develops  completely  before  the  flow  leaves  the  combustion 
chamber,  or  the  flame  front  is  aimed  to  be  far  away  from  chamber's 
walls.  Therefore,  the  increase  of  the  flame  height  necessitates  the 
use  of  larger  chambers.  For  this  reason,  we  consider  the  second 
constraint  as  the  limiting  value  of  the  flame  height. 


z  =  We  (zc  -  Zc,2)  +  (Zh  ~  Zh’A  ~m(ZQ-  Z02 )  ,2g, 

Wc  (ZC,  1  -  zc,2)  +  2MJ  (ZH.l  -  ZH,2)  “  (zO,l  “  zO,l) 

where  z\  is  the  atomic  mass  fraction  of  atom  i  and  the  subscripts  1 
and  2  refer  to  the  value  of  z*  at  the  fuel-  and  air-inlet  streams, 
respectively.  The  flame  height  is  calculated  as  the  maximum  dis¬ 
tance  of  the  iso-value  Z  =  Zstoc  =  0.055  from  the  center  of  the  inlet- 
fuel  nozzle. 


Fig.  18.  The  entropy  generation  rate  due  to  different  phenomena  versus  the  swirl 
number.  The  other  parameters  are  Qf  =  10  1pm,  £  =  2,  Aa  =2  cm2,  </>  =  0.5, 
Va  =  15.876  m/s,  and  Tiniet  =  298.15  K.  Line  legend:  the  same  as  Fig.  9. 
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Fig.  20.  The  entropy  generation  rate  due  to  different  phenomena  versus  the  air-inlet 
velocity.  The  other  parameters  are  Qf  =  10  1pm,  £  =  2,  Aa  =  2  cm2,  S  =  0,  0  =  0.5, 
and  Tiniet  =  298.15  K.  Line  legend:  the  same  as  Fig.  9. 
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The  calculated  flame  height  for  all  cases  of  Figs.  7, 11, 17  and  19  is 
reported  in  Fig.  21.  As  it  is  observed  in  this  figure,  the  bluff-size  has 
the  weakest  and  the  air-inlet  velocity,  especially  for  velocities 
smaller  than  10  m/s,  has  the  strongest  effect  on  the  variation  of  the 
flame  height.  Using  the  results  reported  in  this  figure,  one  can 
choose  a  limiting  value  for  the  flame  height  based  on  their  specific 
size  considerations.  In  our  optimizations,  we  considered  the  value 
of  0.25  m  as  our  second  constraint.  Therefore,  the  optimal  design 


corresponds  to  the  case  with  minimum  entropy  generation  among 
all  cases  which  satisfy  the  following  two  constraints. 

Tmax  E  MAT  —  1900  K 

Flame  height  <  0.25  m  ^ 


5.  Conclusion 


a) 


Va  (m/s) 


In  this  study,  we  simulate  the  turbulent  combustion  of  a  mixed 
bluff-body-swirl  stabilized  flame  in  a  gas  turbine  combustion 
chamber.  The  governing  equations  are  the  continuity,  momentum, 
species  transport,  and  ideal  gas  law.  The  turbulence  closure  is  pro¬ 
vided  by  the  RNG  k-e  model  in  conjunction  with  the  eddy-dissipation 
model  for  the  averaged  reaction  rate.  The  governing  equations  are 
solved  numerically  using  a  finite  volume  method.  Then,  we  use  two 
methods,  i.e.  the  direct  and  indirect  methods,  for  calculating  the 
entropy  generation.  After  grid  independency  check  and  verifying  our 
entropy  generation  computations  by  a  simple  analytical  case,  we 
investigate  the  effects  of  different  parameters,  including  the  equiv¬ 
alence  ratio,  fuel-inlet  flow  rate,  bluff-size  ratio,  swirl  number,  and 
air-inlet  velocity,  on  the  entropy  generation.  We  perform  the  design 
process  of  the  combustion  chamber  by  proposing  the  optimum  value 
of  each  parameter  based  on  the  EGM  method,  i.e.  by  minimizing  the 
total  entropy  generation  under  the  two  maximum  allowable  tem¬ 
perature  and  size  constraints.  We  also  conclude  that  the  two  direct 
and  indirect  methods,  involving  different  approximations,  differ  by 
about  6%  in  average.  The  entropy  generation  rate  due  to  chemical 
reaction  composes  the  largest  portion  of  the  total  entropy  genera¬ 
tion.  The  rates  due  to  heat  transfer,  coupling  between  heat  and  mass 
transfer,  mass  transfer,  and  viscous  friction  have  the  next  places, 
respectively.  And  finally,  the  entropy  generations  due  to  chemical 
reaction  and  heat  transfer  are  opposing  factors  whose  trade-off  re¬ 
sults  in  the  optimum  design  condition. 


Appendix  A 


The  entropy  generation  rates,  Eqs.  (21a-e),  in  axisymmetric 
coordinates  are  written  as 
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Fig.  21.  The  flame  height  for  cases  in  a)  Fig.  7  (Of  =  10  1pm),  b)  Fig.  11,  c)  Fig.  17,  and  d) 

Fig.  19. 


where,  x,  r,  and  0  are  the  axial,  radial,  and  tangential  directions, 
respectively.  Note  that  us  is  called  swirl  velocity. 
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